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A High-Gain Nonlinear Observer 
with Limited Gain Power 

Daniele Astolfi and Lorenzo Marconi 


Abstract —In this note we deal with a new observer for 
nonlinear systems of dimension n in canonical observability form. 
We follow the standard high-gain paradigm, but instead of having 
an observer of dimension n with a gain that grows up to power 
n, we design an observer of dimension 2n — 2 with a gain that 
grows up only to power 2. 

Index Terms —Observability, nonlinear observers, high-gain 
observers. 


I. Introduction 


In this note we consider the problem of state observation 
for nonlinear systems of the form 

z = fiz)+d{t), y = h{z) + iy{t) (1) 


where z G Z C R" is the state, y G R is the measured output, 
/(•) and h{-) are sufficiently smooth functions, d{t) G R” is 
a bounded disturbance and z/ G R is the measurement noise. 
Among the different techniques for observer design available 
in literature (see ifT^ . 16]) we are particularly interested to the 
so-called high-gain methods that have been shown to be effec¬ 
tive in many control scenarios. In this respect we assume that 
the pair (/(•), h(-)) fulfils an uniform observability assumption 
(see Definition 1.2 in ifT^ ). which implies the existence of a 
diffeomorphism f ■. Z ^ R" such that the dynamic of the 
new state variable x = (^{z) is described by the canonical 
observability form (see Theorem 4.1 in 112]) 


X = AnX + B„ip{x) + d{x, t), y = CnX + v{t) , (2) 

where (^(•) is a locally Lipschitz function, 
d(j){z) 


d{x,f) '■= 


dz 


d{t ), 


Z — (f)iriv{x) 


with finvi') the inverse of ff) (namely t^inv ° (j){z) = z for 
all z G Z), and (A„, Bn, Cn) is a triplet in “prime form” of 
dimension n, that is 


An — 
Cn = 


0(n—l)xl In — 1 

0 Olx(n-l) 

1 0ix(„_l) ) . 


Bn = 


0 (, 


.-l)xl 

1 


(3) 


System (|2] is defined on a set A” := (j){Z) C R". 

For the class of systems 0, it is a well-known fact (Q) 
that the problem of asymptotically, in case d{x, t) = 0 and 
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i'{t) = 0, estimating the state x can be addressed by means 
of a high-gain nonlinear observer of the form 

X = AnX + Bn^Psix) + Dn{t)Kn {CnX - y) (4) 

with 

Dn{()=d\ag(^i, , Kn = (^Ci ■ ■ ■ Cn^ , 

where f is a high-gain design parameter taken sufficiently large 
(i.e. £ > i* with £* > 1), the cfs are chosen so that the matrix 
{A -\- KC) is Hurwitz (i.e. all its eigenvalues are on the left- 
half complex plane), and ips{-) ■ K" —?► R is an appropriate 
’’saturated” version of (p{-). As a matter of fact, it can be 
proved that if d{x,t) = 0 and v{f) = 0, if p{-) is uniformly 
Lipschitz in X, namely there exists a (p > 0 such that 

\\p{x')-p{x'')\\<p\\x'-x''\\ Vx',x"gA’, (5) 

and Ps{-) is chosen bounded and to agree with pf) on X, the 
observation error e{t) = x{t) — x{t) originating from (|2] and 
(|4] exponentially converges to the origin with an exponential 
decay rate of the form 

\\x{t) — x(f)|| < exp(—/3ff)||a:(0) — x(0)|| , 

where a and /3 are positive constants, for all possible initial 
condition x(0) G R" as long as x{t) G X. In particular, note 
that the exponential decay rate may be arbitrarily assigned by 
the value of £ with a polynomial ’’peaking” in £ of order n — 1. 
It is worth noting that the uniform Lipschitz condition Q is 
automatically fulfilled if A is a compact set. In case d{x, t) 
or v{f) are not identically zero, as long as they are bounded 
for all t > 0 and for alu x & X, the observer (|4] guarantees 
a bound on the estimation error that depends on the bound of 
d{-, •), of vf) and on the value of £. In particular, the following 
asymptotic bounds can be proved 

lim sup ||x(f) — a;(f)|| < 

t—foo 

7 maxi lim sup ||--r((') d(a;(f), f)|| , lim sup iz(f)|| 1 

I ^ t—^oo \ 

where 7 is a positive constant and 

r(^) = diag(r-\ 1). (6) 

As above the previous asymptotic bound holds for all possible 
i(0) G R" as long as x{t) G X. Note that a high value of 
£ leads to an arbitrarily small asymptotic gain on the n-th 
component of the disturbance d{-,-). On the other hand, a 

'boundedness of d{x,t) is automatically guaranteed if X is compact. This, 
in turn, is the typical case when such observers are used in semiglobal output 
feedback stabilisation problems, tm. 
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large value of ^ is, in general, detrimental for the sensitivity 
of the asymptotic estimate to the sensor noise and to the first 
n — 2 disturbance components. 

Observers of the form © are routinely used in many 
observation and control problems. For instance, the feature 
of having an exponential decay rate and an asymptotic bound 
on the last component of d that can be arbitrarily imposed 
by the value of ^ is the main reason why the above observer 
plays a fundamental role in output feedback stabilisation and 
in setting up semiglobal nonlinear separation principles ( iflTl . 
PI). In that case the set X is an arbitrarily large compact set 
which is made invariant by the design of the state feedback 
stabilisation law and of the high-gain observer. We observe 
that, although the asymptotic gain with respect to v increase 
with the observer is anyway able to guarantee ISS with 

respect to the sensor noise ( lfT4l ). 

The main drawback of observers of the form though, 
is related to the increasing power (up to the order n) of the 
high-gain parameter which makes the practical numerical 
implementation an hard task when n or f are very large. 
Motivated by these considerations, in this note we propose 
a new observer for the class of systems © that preserves 
the same high-gain features of @ but which substantially 
overtakes the implementation problems due to the high-gain 
powered up to the order n. Specifically, we present a high- 
gain observer structure with a gain which grows only up to 
power 2 (instead of n), at the price of having the observer 
state dimension 2n — 2 instead of n. 


II. Main Result 


We start by presenting a technical lemma instrumental to 
the proof of the main result presented in Proposition [T] Let 
Ei G R2x2^ q. g ]g2x2^ i = l,...,n-l,and N € R2x2 
matrices defined as 


E^ = 


-hi 

-ka 



0 hi 
0 ka 



where {hi^ka) are positive coefficients, and let M G 
]^( 2 n- 2 )x( 2 n- 2 ) block-tridiagonal matrix defined as 

/£;iiV0 . 0\ 

Q2 E2 N ■■. : 

0 ■■ ■■■ : 

; ■■■ Q^ E, N \ 


: '. '. '. '. 0 

: ■ ■ Qn-2 En-2 N 

\0 . 0 Qn-1 En-l / 

(7) 

It turns out that the eigenvalues of M can be arbitrarily 
assigned by appropriately choosing the coefficients {hi,h 2 ), 
i = l,...,n — 1, as claimed in the next lemma. 


Lemma 1 Let 'P(A) = -I- miA^”“^ -I-... -I- m 2 „_ 3 A -I- 

?Ti 2 n -2 be an arbitrary Hurwitz polynomial. There exists a 
choice of {hi, ka), * = 1,..., n — 1, such that the character¬ 
istic polynomial of M coincides with V{\). 


The proof of this Lemma is deferred to the appendix where a 
constructive procedure for designing {hi,h 2 ) given V{X) is 
presented. 

The structure of the proposed observer has the following form 

ii = N ^i+iD 2 {£) Kid i = l,...,n-2 

Cn —1 — Cn—l “b E Lps {x ) -j- E 2 {£) K(ji_i'^ Cn—i 

(8) 

where {A, B, C) is a triplet in prime form of dimension 2, 
h G R2, K, = {hi kaf, D2{£) = dmg{£, £^), 

x' = Li^ Li =blkdiag (C, ..., C,/ 2 ) 

-- . -^ 

in- 2 ) times 

C=COl(a,...,en-l) 

ei=y-C^i, d = B'^^i^i-C^i i = 2,...,n-l, 

and (ps{') is an appropriate saturated version of (p{-). 

The variable x' represents an asymptotic estimate of the 
state X of ©. It is obtained by “extracting” n components 
from the state ^ according to the matrix Li defined above. As 
clarified next, the redundancy of the observer can be used to 
extract from ^ an extra state estimation that is 

x" = L 2 ( L 2 = blkdiag {I 2 , B^, B^). 

" -V-^ 

(ra- 2 ) times 

The following proposition shows that the observer ([8]l recovers 
the same asymptotic properties for the two estimates x' and 
x" of the “standard” high-gain observer ©. In the statement 
of the proposition we let 

X = col(a;', x"), x = col(a::, x) . 

Proposition 1 Consider system © and the obsen’er © with 
the coefficients {hi ka) fixed so that the matrix M defined 
in © is Hurwitz (see Lemma\J}. Let Ps{') be any bounded 
function that agrees with ip{-) on X, and assume that d{x,t) 
is bounded for all x C X and for all f > 0. Then there exist 
Ci > 0, t = 1,..., 4, and £* >1 such that for any £> £* and 
for any ^(0) G the following bound holds 

||x(f) - x(t)|| < maxjci exp(-C 2 ^f) ||x(0) - x(0)|| , 

C3\\]md{-)\\oo, C4r-l||K-)l|oo} (9) 

where r(£) is as in ©, for all t > 0 such that x{t) G X. 

Proof: 

Consider the change of variables 
Ci h ■= 6 - col {xi, Xi+i) i = 1,... ,n-l , 
by which system ® transforms as 

it = Hili+N^2-di{t)+D2{£)Kiv{f) 

i = H,l, + Nl^+i+D2{£)K,B'^l^-i-h{t) 

i = 2,... ,n - 2 

in-l = Hn-l in-1 + D2{£) Kn-l B"^ ^n-2 + 

BA(p{i,x) - d„-i{t) 
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where 

Hi = A-D 2 {i) KiC, di(t) = {di{x{t),t),di+i{x{t),t))^ , 

for i = 1 ,..., n — 1 , with di{t) the i-th element of the vector 
d{t), | = col(^i,...,|„_i), and 

x) = + x)- Lp[x) . 

Rescale now the variables as follows 

By letting e = col(ei ... £n-i), an easy calculation shows 
that 

e = iMe + £-(”-1) (B 2 «- 2 A(/ 3 £(e, x) + ve{t) + ne{t)) , 

( 10 ) 

where B 2 n -2 is the zero column vector of dimension 2n — 2 
with a 1 in the last position, and 

viit) = -col {rD 2 ii)-^di, r-^+^D 2 ii)-^d,, 

e^D2{i)-^dn-l) , 

= i'^Kiv{t) , Ki = col(ifi, 0,..., 0) , 

and 

A(fie{e,x) = (ps{LiS{£)e + x) - ip(x) , 

where S{£) = diag(i , 1, £, ..., (g) D 2 {£). Being ip{-) 

uniformly Lipschitz in X and (ps{-) bounded, there exists a 
(5i > 0 , (52 > 0 and (53 > 0 such that 

\\£-(-^-^) A(pi>{£,x)\\ < (5i||e||, 

< d2\\£Dni£)-^d{xit),t)h 
\\£-^-^^n,{t)\\ < h\\£m\. 

for all e: S x & X and £ > 1. The rest of the 

proof follows standard Lyapunov arguments that, for sake of 
completeness, are briefly recalled. Let P = be such that 
PM + Xl'^P = —I and consider the Lyapunov function 

V = Pe. Taking derivative of V along the solutions of 
(fTOl i. using the previous bounds and letting £* = 4(5i||P||, one 
obtains that there exist positive constants oi, 02 , 03 such that 
for any £> £* 

||e|| > max{ai||i:)n(f)“^d(a;(f),f)||,a 2 ||i^(f)||} =?> 

V <-a:i£\\ef. 

As P is symmetric and definite positive, it turns out that 
Aik IP A y{s) A ^IklP where A and A are respectively 
the smallest and the highest eigenvalue of P. By using these 
bounds, the previous implication leads to conclude that, as 
long as 

V >X max {aj \\Dn{£)~^ d{x{t), f)|p , al |k(f)||} 

then V{t) < exp(—^f)T4(0). By using again the bound on 

V in terms of A and A, the following estimate on e(f) can be 
easily obtained 

||e(f)|| < max{ a 4 exp(-a 5 £t) ||e( 0 )|| , 

06 \\Dn{£)~^ d{x{t),t)\\ , a7|k(f)||} 


where 04 = \/X/X, 05 = az/2.X, 05 = \/A/Aai, 
07 = \/X/Xa 2 - Now, using the fact that, for all > 1, 
£-{n-i) Ill'll ^ ii^ii £ iilii^ j-jjg previous bound leads to the 
following estimate on ^(f) 

|||(f)|| < max{a4r“^exp(-a5ff) |||(0)|| , 

06^"^ \\Dn{£)~~^ d{x{t),t)\\ , a7^”"kk(f)||} 

by which the claim of the proposition immediately follows by 
bearing in mind the definition of r(£), x, x, and by noting 
that Hill < ||x-x|| < 2|||||. ■ 

III. About the sensitivity of the observer to high 

FREQUENCY NOISE IN THE LINEAR CASE 

The trade-off between the speed of the state estimation and 
the sensitivity to measurement noise is a well-known fact in the 
observer theory. In this respect, high-gain observers tuned to 
obtain fast estimation dynamics are necessarily very sensitive 
to high-frequency noise. Bounds on the estimation error in 
presence of measurement noise for the standard high-gain 
observers have been studied, for instance, in ifTSl and Q, and 
different techniques have been developed in order to improve 
rejection, mainly based on gain adaptation (see, among others, 

El, El). 

In this section we compare the properties of the standard 
high-gain observer (HI and the proposed observer ® with 
respect to high-frequency measurement noise by specialising 
the analysis to linear systems. 

In particular we consider systems of the form ^ with (p{-) 
a linear function of the form (p{x) = $ x where $ is a row 
vector of dimension n. Moreover, in this contest, we consider 
d{t) = 0 and 

iy{t) = aNsiii{ujNt +Jn) (H) 

where qn > 0, > 0, and /at are constants. It is shown 

that the ratio between the asymptotic estimation error on 
the i-th state variable provided by the new observer (O and 
the one provided by the standard observer (HI is a strictly 
decreasing polynomial function of the noise frequency for 
i = 2,...,n. In this regard the new observer has better 
asymptotic properties with respect to high-frequency noise as 
far as the state estimation variables are concerned (except for 
the first one). This is formalised in the next proposition. 

Proposition 2 Let n > 1. Let Kn G M”, Ki G i = 
1,..., n — 1 and £ G R fixed so that the error dynamics 
of the observers © and are Hurwitz- Moreover, with p G 
{1, 2,..., n} denoting the position index of the first non zero 
coefficient of the vector $ (and p = n if ^ is the zero vector), 
let r[ > 1 be the constants defined as 

r- = min {i , (n — 1) , {p + n — i + l)} i = 1,... ,n . 

There exist ujfi > 0 and Ci > 0 such that for all lon > 

On > 0 and /at we have 

limsup |x'(f) — Xi{t)\ 

^ _ —(r' —1) 

limsup \xi{t) — Xi{t)\ ~ 


V * = 1,..., n. 



4 


Proof: Consider system © and the standard high-gain 
observer (|4|. By letting e = £ (x—x), the e-dynamics 
read as 


e = £(A^ + K„C„) e + 0„(£) e + £K„iy(t) (12) 

where 0„(£) = —£>„(£) with D„(£) defined in (|4]i. It is a 

£n 

linear system that is Hurwitz by the choices of Kn and £. 
Similarly, consider system (|2]i and the new observer (| 8 ]l. With 
£ defined as in the proof of Proposition 1 we obtain system 
(doll compactly rewritten as 

£ = £M£ + B2n-2$ On(£)Li £ + £Kiv{f) . (13) 


It is an Hurwitz system by the choices of Ki, i = 1,..., n — 
1, and £. We consider now the n systems given by the 
dynamics (fT 2 li with input v and with outputs 

Xi — Xi = t~^ei, i = 1,.. .n , 

and we denote by : R —^ C, i = 1,..., n, the harmonic 
transfer functions of these systems. A simple computation 
shows that these systems have relative degree = 1 , for all 
i = 1,... ,n. Similarly, we consider the n systems given by 
the dynamics (fTsT l with input z/ and outputs 

x^ — Xi = P~^£i^i i = l,...n — l, 

^ ^n— 1,2 ■) 

and we denote by F/ : R —>■ C, i = 1,..., n, the harmonic 
transfer functions of these systems. Simple computations show 
that these systems have relative degree r', for alH = 1 ,..., n, 
with r' defined in the statement of the proposition. By defini¬ 
tion of harmonic transfer function and by the fact that systems 
(fT^ and (foi l are Hurwitz, it turns out that 

limsup - a;i(f)| = \Fi{u}N)\aN 

limsup = |F/(wAr)|aAr 


for any ujn > 0 and Vi = l,...,n. Furthermore, by the 
fact that Fi{uj) and Fl{uj) have, respectively, relative degrees 
ri = 1 and r[, i = ..., n, it turns out that there exist positive 
Ci, c' and ujfj > 0 such that 


|Fi(w)| > CiOJ ^ 

|F/(c^)| < 


yu}>uj*^, 


by which the result immediately follows. 


IV. Example: Observer eor the uncertain Van Der 
Pol Oscillator 

Let consider the uncertain Van der Pol oscillator 

z = —a^z -I- /3(1 — , y = z + v{t) , (14) 

where y € R is the measured output and a, /3 are uncertain 
constant parameters. We let y = {a^, /3)^ and we assume that 
y & U, with U a compact set of R^ not containing the origin. 
The state (z, z) belongs to a compact invariant set W C R^, 
which is the limit cycle of the Van Der Pol oscillator. We 
observe that W depends on p. Following 0 , system (Ell 
extended with fi = 0 can be immersed into a system in the 
canonical observability form (l2]i with n = 5. As a matter of 


fact, let := col {z^'‘\ ..., z^A'j be the vector of time 

derivatives of z, with 0 < i < j, and let X be the compact set 
of R^ such that Z[q i] G W Z[o, 4 ] G X. Simple computations 
show that 

^[2,4] = Tf'(^[0.3]) F , = P(^[0.4]) P 

where 

(-Z {l-z^)z \ 

^(^[o.a]) = -^ z- 2 yp-z‘^'i 

y—z z^^^ — 2 z^ — 6 zzz — z^z^^^ j 

and 

p(' 2 [o, 4 ]) = {z^^^ , z^‘^\l — z^) — 12 z^z — 6 zz^ — 8 zzz(^^) 

with (see 0 ) 

rankT(z[o, 3 ]) = 2 for all Z[o, 4 ] G X . 

Hence, by letting 

¥’(^[ 0 . 4 ]) = p{z[ 0 ,i]) Kz[0,3]) 

with 

Kz[0A]) = ^^'*'(^[ 0 . 4 ]) 2 :[ 2 , 4 ] , 

where T’l'(-) is the left-inverse of T(-), and letting x = -Z[o, 4 ], 
it is immediately seen that system (fT4li and fi = 0 restricted 
to W X W is immersed into the system 

X = A 5 X + B 5 (p{x), y = Cs x + r/(f) , (15) 

where is a triplet in prime form of dimension 5. 

By following the prescriptions of Section II, we imple¬ 
mented the proposed observer ([^ as 

a = A^i+N^2 + D2{£i)K^{y-C^i) 

= A^i + By:s{.x') + D2{,£i)Ki{B^^^-CU) 

(16) 

where ips{-) is any locally Lipschitz bounded function that 
agrees with (/?(•) on X, = (^a,Cz 2 ) and the coefficients 
of Ki are ku = 0.6, fci 2 = 0.3, ^21 = 0.6, ^22 = 0.111, 
fcai = 0.6, ^32 = 0.04 85, k 4 i = 0.6, £142 = 0.0178, such 
that the roots of V{X) are —0.1, — 0 . 2 , — 0 . 2 , —0.3, —0.3, 
—0.4, —0.4, —0.5. With the same notation of (01 we have 
C = ( 6 , 6 , 6 ,^ 4 ), x' = Li^ and f" = L 2 C 

In the simulations we fixed a = 1, /? = 0.5, gain £1 = 100 
and initial conditions (z,z) = (1,0) for (fT4l i and = 0 for 
(Ell- Figures 0 and 0 show the error state estimate \x' — x\ 
and \x" — x\ of the proposed observer (El for the first two 
components (namely the estimation of the state of (El) when 
there is not sensor noise. 

By following Section III we compared the observer 
(El with a standard high-gain observer in presence of 
high-frequency sensor noise, numerically taken as v{t) = 
10“^ sin(10^ f). The high-gain observer has been implemented 
as 

X = A 5 X + B^ifsix) + D{£2)K^(%j - C5X) (17) 

where K 5 = (1.5, 0.85, 0.225, 0.0274, 0.0012)'^ so that 
the eigenvalues of (A 5 + K^Cf) are —0.1, —0.2, —0.3, 
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—0.4, —0.5, and £2 = 100. Table 1 shows the normalized 
asymptotic error magnitudes of the proposed observer 
and the standard high-gain observer (E), where the 
normalized asymptotic error for the i-th estimate is defined as 

11 — ati 11 a = lim sup \\xi{t) — Xi{t)\\/aN- Although the result 

t—^OO 

in Section III is given just for linear systems, the numerical 
results shown in the table show a remarkable improvement 
of the sensitivity to high-frequency measurement noise of the 
new observer with respect to the standard high-gain observer. 


Standard High Gain 

Modified 


Modified 

Observer £ 

Observer x' = Li^ 

Observer £ 

" = ^2? 

||£i - xi\ 

a =0.15 

ll^j - XlWa = 0.06 

1 

*j' - 1 

a = 0.06 

*2 - X2 \ 

a =8 

11^2 “ = 0.2 

1 

£ii — xi 1 

a = 3 

(ra - 

a = 2 ■ 102 

11*3 - ^aiia = 0.2 

1 

£3 — X3 1 

a = 3 

£4 - xf 

a =2.5- 103 

11*4 - a:4iia = 0.1 

1 

£4 - xf 

a = 2 

(rs - 2:5 

a = 104 

il*5 - *5i|a = 0.3 

1 

*5 -*5 

a = 0.3 


Table 1: Normalized asymptotic errors in presence of noise. 


numerically simulated in the previous section. A complete 
characterisation of the sensitivity to sensor noise of the new 
observer is an interesting research topic are that is now under 
investigation. 

The peaking phenomenon due to wrong initial conditions 
and fast convergence that is typical of high-gain observers 
is not prevented by our proposed structure. However, other 
techniques to deal with peaking (such as saturations, time- 
varying gains M, gradients techniques El, and others) 
are available and can be adopted to improve the proposed 
observer structure. In this work we didn’t consider the multi¬ 
output case. For the specific class of multi-output systems 
which are diffeomorphic to a block triangular form in which 
each block is associated to each output and it has a triangular 
dependence on the states of that subsystem (see El), the 
proposed structure can be simply applied block-wise to obtain 
a high-gain observer. Apart this case, a complete extension to 
the multi-output case is not immediate and under investigation. 




Fig. 2; Error state estimate \x 2 — X 2 \ (blue line) and — X 2 \ 
(red line). 


V. Conclusions 

We presented a new observer design based on high-gain 
techniques with a tunable state-estimate convergence speed. 
With respect to standard high-gain observers the state dimen¬ 
sion is larger (2n — 2 instead of n) with a clear benefit in 
the observer implementation due to the power of the high- 
gain which is only 2 and not n. Moreover, when specialised 
to linear systems, we showed the benefit of the proposed 
observer with respect to the standard high-gain observer in 
terms of high-frequency noise rejection. Benefits that are 
clearly confirmed also for the nonlinear Van-der Pol example 


Acknowledgment. We wish to thank Laurent Praly for 
suggesting the design procedure presented in Appendix A. 

Appendix 

A. Procedure to assign the eigenvalues of M 

Consider the matrices Mi € x recursively 

defined as 

Mi=Ei, M, = 2)’ * = 2,...,n-l 

where iV* = col(02(j_2)x2, N), Qi = {02{i-2)x2, Qi) and Ei, 
i = 1,... ,n — 1, Qi, i = 2,... ,n — 1, and N are defined as 
in the definition of M. Note that M = Mn-i and, by letting 
Ki = {kii ki 2 )'^, note that Qi and Ei depend on Ki, while Mi 
depends on Ki,... ,Ki. We let Vm^^) = -|- 

... -I- rn\^_^\ + and 7^Mi_i (A) = A^*“^ -b -b 

.. •+TO 27 - 3 A-bm 27.^2 characteristic polynomials of Mi and 
Mi-i, and we use the notation = col(mj,..., rrij) G 

RC = co\{m]~^,..., G R^ for some j < 2i and 

k <2i-2. 

The characteristic polynomial of VMi (A) is computed as 
^Mi(A) = A(A-b (A) 

+ki2 [7^Mi_i (A) — A(A-b (A)] . 

Hence, simple, although lengthy, computations show that 
the coefficients 2 i] of TmiW and TO[ 72 i- 2 ] of 'PMi_i(A) 
are related as follow 

"^[1.2i-2] = (-f2*-2 + fcilF')TOji"2i-2] + 

TO2*-1 = ^*1^2772 (18) 

mli = h2 

where vi G R^®“^ is the zero vector with a 1 in the first 
position, and E G R( 2 *- 2 )x( 2 z- 2 ) ^ero matrix with the 

identity matrix Qi-s in the lower left block. Note that (/ 2 i -2 + 
kiiE) is invertible for all kn. Hence, from the first of (fTsT l. 
one obtains 

'”^[1,21-2] “ ^(^[l,2i-2]i 
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where 

—2] 1 ^il) “ (-^2i—2 + knF) (W[1^2i-2] ~ Vi ) , 

which, embedded in the second and in the third of (fTST l, yield 
the relations 


where 

Cr2(TO[i^2i]’ ^il) = 


0, ki2 = a2{rn\^ 2i]^ kii) 

hi vj A(m|^_2i-2]’ “ ^2*-l 

»^2z 

vl’ A(m^^ 2._2], hi) 


in which V 2 G is the zero vector with a 1 in the last 

position. We observe that (Ti(-, •) is a polynomial in hi of odd 
order 2i — 1. As a consequence, for any 2 i_i] there always 
exists at least one real hi fulfilling ai{m)-^ 2 z-i]’ ^*i) “ 

The previous results can be used to set up a ’’basic as¬ 
signment algorithm” that is then used iteratively to solve the 
eigenvalues assignment of the matrix M. 

Basic assignment algorithm. Let 'Pi(A) = + 

... + m 2 i be an arbitrary polynomial. Then, there exist a 
real Ki = {hi,ki 2 )'^ and a polynomial Vi-i{X) = A^®“^ + 
+ ... + fh) 2 il 2 ®*tch that 


K, = K, 
'PMi-i = 'Pi-1 (A) 


Pm,(A) = P.(A). 


As a matter of fact, by letting 2i_i] the coefficients 
of Pi{X), it is possible to take hi as a real solution of 
^l(Ai[1.2i-l]’^il) = 0, fci2 = Cr2(fhJ;^ 2i]> ^*l)’ ™d to 
take the coefficients ^ 1 ^^ 121 - 2 ] the polynomial Pi_i(A) as 
™[l, 2 i- 2 ] “ A(frijj^ 2 i- 2 ] ’ ^*1) ■ ^ 

With the previous algorithm in hand, the design of 
Ki, ..., Kn-i to assign an arbitrary characteristic polynomial 
to M, can be then immediately done by the following steps; 

1) With 'P„_i(A) the desired characteristic polynomial of 
M, compute {Kn-i,Vn- 2 {X)) by running the basic 
assignment algorithm with i = n — 1. 

2) Compute iteratively {Ki,Pi-i{X)) by running the basic 
assignment algorithm for i = n — 2,... ,2. 

3) Compute Ki = iku.ka)'^ so that A^ + fc^A + ka = 
Pi (A). 
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